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Abstract. Evidence is presented from numerical magneto-hydrodynamical simulations for the existence of magnetic 
activity in late-type giant stars. A red supergiant with stellar parameters similar to that of Betelgeuse {a Orionis) 
is modeled as a "star-in-a-box" with the high-order "Pencil Code" . Both linear kinematic and non-linear saturated 
dynamo action are found: the non-linear magnetic field saturates at a super-equipartition value, while in the linear 
regime two different modes of dynamo action are found. It is speculated that magnetic activity of late-type giants 
may influence dust and wind formation and possibly lead to the heating of the outer atmospheres of these stars. 
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1. Introduction 

There are indications from both dynamo theory and ob- 
servations that some late-type giant stars such as red su- 
pergiants and asymptotic-giant-branch stars (AGB stars) 
may harbor magnetic fields. On the theoretical side, it has 
been suggested that non-spherically symmetric planetary 
nebulae (PNe) may be a result of the coUimating effect of 
a strong magnetic field: Blackman et al. 1)200 l|l 



studied in- 
terface dynamo models similar to the mean field theory's 
solar aw-dynamo and found that the generated magnetic 
fields were strong enough to shape bipolar outflows, pro- 
ducing bipolar PNe, while also braking the stellar core 
thereby explaining the slow rotation of many white dwarf 
stars. Also using mean field dynamo theory Soker & Zoabi 
H2002|l propose instead an a'^uj dynamo due to the slow 
rotation of AGB stars rendering the w-effect ineffective. 
They find that the magnetic field may reach strengths 
of ~ 100 Gauss, significantly less than that found by 
Blackman et al. (|2001|l . On the one hand, they believe 
that the large-scale field is strong enough for the forma- 
tion of magnetic cool spots (see also Soker & Kastner ,200BI 
on AGB star fiaring). These spots in turn may regulate 
dust formation, and hence the mass-loss rate, but the au- 
thors argue that they cannot explain the formation of non- 
spherical PNe (see also Soker|2)02): on the other hand, the 
locally strong magnetic tension could enforce a coherent 
flow that may favor a maser process. 



On the observational side of things, maser polarization 
is known to exist in circumstellar envelopes of AGB stars 
(e.g. Gray et al. 119991 Vlemmings et al. 120031 and recently 
Sivagnanam 1200 j|l and X-ray emission has been observed 
from some cool giant stars (e.g. Hiinsch et al. 119981 and 
Ayres et al. l2008|l . These observations are generally taken 
as evidence for the existence of magnetic activity in late- 
type giant stars (cf. Soker & Kastner 2003). 

The cool star Betelgeuse (a.k.a. a Orionis) is an ex- 
ample of an abundantly observed late-type supergiant 
that displays irregular brightness variations interpreted 
as large-scale surface structures (e.g. Lim et al. ^98 and 
Gray 2000 ) . It is one of the stars with the largest apparent 
sizes on the sky — corresponding to a radius in the interval 
600-800 Rq. Freytag et al. (20(2|) performed detailed nu- 
merical 3-d radiation-hydrodynamic (RHD) simulations of 
the convective envelope of the star under realistic physical 
assumptions, while trying to determine if the star's known 
brightness fluctuations may be understood as convective 
motions within the star's atmosphere: the resulting mod- 
els were largely successful in explaining the observations 
as a consequence of giant-cell convection on the stellar sur- 
face, very dissimilar to solar convection. Dorch & Freytag 
(2002) performed a kinematic dynamo analysis of the con- 
vective motions (i.e. not including the back-reaction of the 
Lorentz force on the flow) and found that a weak seed mag- 
netic fleld could indeed be exponentially amplified by the 
giant-cell convection on a time-scale of about 25 years. 



This paper reports on full non-linear magneto- 
Send offprint requests to: S.B.F. Dorch — dorch@astro.ku.dk hydrodynamical (MHD) numerical simulations of dynamo 
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action in a late-type supergiant star with fundamental 
stellar parameters set equal to that of Betelgeuse. The 
paper is organized as follows: Section[21contains a descrip- 
tion of the numerical model, code and setup, Section|3|de- 
scribes and discusses the results, the convective flows and 
dynamo action, and Section 0] contains a short summary 
and conclusion. 



2. Model 

The full 3-d MHD equations are solved for a fully convec- 
tive star. This is an example of a "star-in- a-box" simula- 
tion, where the entire star is contained within the com- 
putational box. The computer code used is the "Pencil 
Code" by Brandenburg & Dobler-'^ . The code has been em- 
ployed in several astrophysical contexts including e.g. hy- 
dromagnetic turbulence (see Brandenburg & Dobler 2002T 
Dobler et al. 120031 and Haugen et al. l20()3|l . The numerical 
method performs well on many different computer archi- 
tectures especially on MPI machines, and uses 6th-order 
spatial derivatives and a 2A^-type 3rd-order Runge-Kutta 
scheme. The code has a "convective star" module that al- 
lows the solution of the non-linear MHD equations by the 
numerical pencil scheme in a star with a fixed radius R 
and mass M. The code solves the following general form 
of the compressible MHD equations: 



Ding 

m 

Du 

Dt 



OA 

'dt 

,Ds 

m 



- -V u, 

Q_ / s \ i X B 

= -civ — + In - V$ + 

Vcp / g 

+ V ^V^u + i VV ■ u + 25 ■ V In £1 
= u X B - 7?//oj, 

= H-C + V ■ (XVT) + r^^iof + 2^?^.S^ 



(1) 

(2) 

(3) 
(4) 



where Eq. (1) is the mass continuity equation and Eq. (2) 
is the equation of motion. Here g is density, u the fluid 
velocity, t is time, D/Dt = d/dt -f u ■ V is the comoving 
derivative, cs is the sound speed, s is the entropy, $ is the 
gravity potential, j = V x B//io the electric current den- 
sity, B the magnetic flux density, v is kinematic viscosity, 
and S is the traceless rate-of-strain tensor. Eq. (3) is the 
induction equation where A is the magnetic vector poten- 
tial, B = V X A the magnetic flux density, rj — 1/ {fioa) is 
the magnetic diffusivity {a being the electrical conductiv- 
ity) , and fiQ the magnetic vacuum permeability. Eq. (4) is 
the energy, or rather entropy equation, where T is tem- 
perature, Cp the specific heat at constant pressure, 7i and 
C are explicit heating and cooling terms, and K is the 
thermal conductivity (radiation is not taken into account 
in MHD). 

Variables are computed in terms of R and M so 
that e.g. the unit of the star's luminosity L becomes 
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^(GM/R)^. In the present case these fundamental pa- 
rameters are set to R = 640 Rq and M = 5 Mq yielding 
a luminosity of L = 46000 L0, consistent with current 
estimates of Betelgeuse's size, mass and luminosity. The 
model employs a fixed gravitational potential <i>, an inner 
tiny heating core entering into Eq. (4) through 7i, and an 
outer thin isothermally cooling spherical surface at r = R 
(corresponding to C in Eq. 4), with a Newtonian cooling 
time scale set to TcooI = 1 year corresponding to the typ- 
ical convective turn-over time in the model of Freytag et 
al. ijMa . 

Due to constraints on computer time and resolution 
the true thermodynamic range of the star can not be rep- 
resented and the surface is cooled at a temperature that 
is 4.5 times higher than Betelgeuse's effective tempera- 
ture which is about Teff = 3500 K (various estimates exist 
in the literature, see e.g. Freytag et al. 2002). The fixed 
gravity is approximately given by a 1/r-potential (as in 
Freytag et al. I2002|l and initially the star expands result- 
ing in a slight decrease of 1.4% of the mean mass density, 
subsequently it re-contracts by roughly 0.1%. Betelgeuse is 
only slowly rotating and a rotational frequency was chosen 
corresponding to a surface rotational velocity of 5 km/s, 
yielding a large Rossby number. 

Models with different numerical resolutions have been 
run for testing: in this paper the results come from a 
model with 128'^ uniformly distributed grid points yielding 
a spatial resolution of Ax = I5R0 the physical size of the 
box being R"^. The models were computed at the Danish 
Center for Scientific Computing Horseshoe 512 cpu Linux 
cluster typically allocating between 16 and 32 cpus. 

Boundary conditions on the computational box are 
anti-symmetric for components of the vector fields u and 
B in the direction of the component (yielding a vanish- 
ing value at the boundary) and symmetric across the 
boundary in the direction perpendicular to the compo- 
nent (yielding a vanishing gradient across the bound- 
ary). Boundary conditions for the density Inp are anti- 
symmetric with respect to an arbitrary value across all 
boundaries (yielding a vanishing second derivative), and 
the boundary condition for the entropy s/cp corresponds 
to a constant temperature at the boundary. 



2.1. Magnetic Reynolds number 

Dynamo action by flows are often studied in the limit of in- 
creasingly large magnetic Reynolds numbers Rcm — iV/r], 
where £ and U are characteristic length and velocity scales. 
Most astrophysical systems are highly conducting yielding 
small magnetic diffusivities i], and their dimensions are 
huge resulting in huge values of RCm- Betelgeuse is not an 
exception and most parts of the star is better conducting 
than the solar photosphere that has a magnetic diffusivity 
of the order of 77 « 10^ m^/s: 

Figure ^ shows the average Spitzer's resistivity as a 
function of radius in the model of Betelgeuse by Freytag et 
al. (,2002,1 . Spitzer's formula (e.g. Schrijver & Zwaan i20()()|l 
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Fig. 1. Diffusivity: average radial magnetic difFusivity rj 
(m^/s) in the model of Betelgeuse as a function of radial 
distance in solar radii Kq . Adopted from Dorch & Freytag 
[MM 

assumes complete ionization and hence the precise values 
of 7] are uncertain in the outer parts of the star, where the 
atmosphere borders on neutral. There is some uncertainty 
connected also with defining the most important length 
scale of the system, but preliminarily taking £ to be 10% 
of the radial distance R from the center (a typical scale of 
the giant cells), and U = urms along the radial direction 
yields Rcm = 10^°-10^^ in the interior part of the star 
where R < 700 Rq. 

In the present case we cannot use that large val- 
ues of Rem (partly due to the fact that we are employ- 
ing a uniform fixed diffusivity 77), but rely on the re- 
sults from generic dynamo simulations indicating that re- 
sults converge already at Reynolds numbers of a few hun- 
dred (e.g. Archontis, Dorch, & Nordlund 12003^ [M)3a^) . 
Furthermore, Dorch & Freytag H2002(l obtained kinematic 
dynamo action in their model of a magnetic Betelgeuse 
at Rcm ~ 500. Here we use an 77 that is about 10® 
times too large compared to the estimated surface value 
in Betelgeuse, leading to a magnetic Reynolds number of 
Rcni ~ 300 (based on the largest scales). 

3. Results and discussion 

There is some disagreement as to what one should require 
for a system to be a "true" astrophysical dynamo. Several 
ingredients seem to be necessary: the flows must stretch, 
twist and fold the magnetic field lines (e.g. Childress & 
Gilbert ri995|l : reconnection must take place to render the 
processes irreversible; weak magnetic field must be cir- 
culated to the locations where flow can do work upon 
it (cf. Dorch l?000|l : and flnally, the total volume mag- 
netic energy Emag = /y CmdV must increase (the linear 
regime) or remain at a constant saturation amplitude on 
a long time scale (non-linear regime). These points are 
based largely on experience from idealized kinematic and 
non-linear MHD dynamo models; e.g. Archontis et al. 
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H2003bl I2003aj) . This paper deals mainly with the ques- 
tion of the exponential growth and saturation of Emag. 




Fig. 2. Simulated surface intensity snapshots at four dif- 
ferent instants, time = 256, 347, 457 and 494 years (from 
upper left to lower right). 



3.1. Convective flows 

Although not the main topic of this paper, it is in its 
place to discuss here also the properties of the convective 
flows in the model, since these supply the kinetic energy 
forming the basic energy reservoir for any dynamo action 
that might be present. It is not expected that the flows 
match exactly what is found in more realistic RHD sim- 
ulations, but at least a qualitative agreement should be 
inferred since the fundamental parameters of this MHD 
model and the RHD model of Freytag et al. I|2002|l are 
the same. 

The velocity is initialized with a random flow with a 
small amplitude. Rapidly large-scale convection cells de- 
velop through-out the star: the giant cell convection is ev- 
ident in both the thermodynamic variables, such as tem- 
perature and gas pressure, as well as in the flow fleld, 
but the observational equivalent however, is the surface 
intensity. Since the model does not incorporate realistic 
radiative transfer (as opposed to the model of Freytag et 
al. l2002|l . only a simulated intensity can be derived. Using 
a frequency independent LTE source function a simulated 
intensity Igim can be defined as: 

Isim(y,z) = rf7T(x,y,z)4e-^('^'y'^)dT(x,y,z), (5) 
Jo 
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t(x) = / Kop(x,y,z)dx, (6) 
Jo 

where T is temperature, r is a measure of optical depth 
along one of the axes (here the radial line of sight direc- 
tion is taken to be the grid x-axis), p is mass density, kq is 
a constant gray opacity and a Stefan's constant. Figure 
12 shows simulated intensity snapshots at four different 
instances; the typical contrast between bright and dark 
patches on the surface is 20-50%, and only 2-4 large cells 
are seen at the stellar disk at any one time corresponding 
to a hand full of cells covering the entire surface. This is 
in qualitative agreement with the RHD models of Freytag 
et al. H2002fl : the surface is not composed of simply bright 
granules and dark intergranular lanes in the solar sense — 
sometimes the pattern is even the reverse of this — e.g. in 
Figure|21the simulated intensity snapshot at time t — 457 
years, the cool dark area in the center of the stellar disk 
is actually a region containing an upward flow. 
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Fig. 3. Powerspectrum of u^(A:) at time = 457 years, and 
a line corresponding to Kolmogorov scaling (dashed line) . 
The vertical lines, from left to right, denote the wavenum- 
bers corresponding to the computational box size fco, the 
stellar radius fcpi, and the Nyquist frequency fcNy (the nu- 
merical resolution). 

More quantitatively, the kinetic power-spectrum (Fig. 
O illustrates that there is much more power on large scales 
than on the small scales of the velocity field: below a 
wavenumber of A; ^ 20 (based on the box size in units 
of the star's radius) power is decreasing fast, but at larger 
scales the power is proportional to k~^^^ corresponding 
to normal Kolmogorov scaling, the inertial range spans 
however only roughly one order of magnitude. In conclu- 
sion the large-scale convective patterns are then typically 
larger than 15-30% of the radius, and are actually often 
on the order of the radius in size. The corresponding ra- 
dial velocities range between 1-10 km/s in both up and 
down flowing regions. 

There are at least three different evolutionary phases 
of convection in the simulations, depending on the level 
of the total kinetic energy Ekin of the convection motions: 



initially there is a transient of about 30 years after which 
the RMS velocity field reaches a level where it fluctuates 
around a value of about 800 m/s (this corresponds to the 
kinematic phase of the dynamo, where the flow is unaf- 
fected by the presence of the still weak magnetic field, see 
below). During the rest of the simulation after about 290 
years, the RMS speed measured in the entire box decreases 
to 500 m/s (when the energy in the magnetic field be- 
comes comparable to the kinetic energy density). During 
the entire simulation, however, the maximum speed in the 
computational box fluctuates around a constant value of 
about 90 km/s. The flows are not particularly helical and 
the mean kinetic helicity is on the order of 10~^ m/s^. 
Mean held acj-type solar dynamos do not produce large- 
scale fields if the kinetic helicity is less than a certain value 
(cf. Maron & Blackman l2002|l and hence we cannot expect 
a large-scale toroidal field in the solar sense to be gener- 
ated. 
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Fig. 4. Linear regime: energy as a function of Betelgeusian 
time in years. The upper almost horizontal line is total 
thermal energy Eth (dashed-dotted), middle curve with 
wiggles is total kinetic energy Ekin (dotted), and lower full 
curve is Emag- The three thin dashed lines corresponding 
to exponential growth with characteristic time-scales of 
3.8, 25 and oo years. 



3.2. Dynamo action 

In an earlier kinematic study of Betelgeuse using a com- 
pletely different numerical approach (Freytag et al. 120021 
and Dorch & Frevtag f2002|l . dynamo action was obtained 
when the specified minimum value of RCm was larger than 
approximately 500 and at lower values of Ke^ the to- 
tal magnetic energy decayed. In the present case RCm 
is of the same order of magnitude and we find an ini- 
tial clear exponential growth over several turn-over times, 
and many orders of magnitude in energy. Figure^ shows 
the evolution of Emag as a function of time, for the first 
225 years (in Betelgeusian time): initially there is a short 
transient, where the field exponentially decays because the 
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Fig. 5. Transition to the non-linear regime: RMS mag- 
netic field B in the whole computational box in Gauss as 
a function of Betelgeusian time in years. The upper curve 
is the equipartition field strength Boq corresponding to the 
average kinetic energy density of the fluid motions (dot- 
ted curve) and the lower full curve is the actual RMS field 
strength B (full curve). The dashed thin curve correspond 
to growth times of 25 years and a horizontal reference line 
at a field strength of 100 Gauss. 



fluid motions has not yet attained their final amplitude. 
Once the giant cell convection has properly begun how- 
ever, the magnetic field is amplified and we enter a linear 
regime of exponential growth. There are two modes of 
amplification in the linear regime; the initial mode with a 
growth rate of about 4 years, which in the end gives way 
to a mode with a smaller growth rate corresponding to 25 
years. This is a slightly unusual situation, since normally 
modes with smaller growth rates are overtaken by modes 
with larger growth rates (cf . Dorch 2000) ; the explanation 
is that while both modes are growing modes, only the one 
with the largest growth rate is a purely kinematic mode — 
while the exponential growth of the second mode is linear 
it is not kinematic — the presence of the magnetic field is 
felt by the fluid through the back-reaction of the Lorentz 
force in Eq. (2) becoming important. This quenches the 
growth slightly and henceforth one can refer to the second 
mode as a "pseudo-linear" mode. 

No exponential growth can go on forever and even- 
tually the magnetic energy amplification must come to a 
halt: the question is then whether the magnetic field re- 
tains a more or less constant saturation value, or if it dis- 
sipates. The latter is only possible if the non-linear mode 
is a decaying mode that corresponds to a negative growth 
rate. In case of saturation the typical field strength is ex- 
pected to be on the order of the equipartition value corre- 
sponding to equal magnetic and kinetic energy densities. 
Figure El shows the RMS magnetic field strength Brms 
within the entire model star as a function of time for ~ 700 
Betelgeusian years: the pseudo-linear mode as well as the 
mode in the non-linear regime are shown. The RMS mag- 
netic field saturates at a value slightly above the RMS 



equipartition field strength Bcq = \J (pu^) ^ 90-100 
Gauss, corresponding to a value of about 120-130 Gauss. 
In terms of total energy this means that the magnetic en- 
ergy E„iag is above equipartition with the kinetic energy 
Ekin by roughly a factor of two. Hence the field cannot 
be said to be extremely strong, but it is not particularly 
weak in most parts of the star either. What may be inter- 
esting from an observational point of view is the strength 
of the field at the surface. In the non-linear regime the 
field strength at the sphere with radius r = R can be up 
to ~ 500 Gauss, while in the interior of the star it can be 
as high as a few kG: the strong field of the intermittent 
magnetic structures almost completely quenches the ve- 
locity field in these regions that are small-scale compared 
to the scale of the convection; i.e. the local field can be far 
above equipartition. 
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Fig. 6. Distribution of magnetic field strength: PDF for 
the magnetic field |B| (in kG) at time = 732 years within 
the non-linear saturated regime. 



3.3. Magnetic structures 

It is interesting to examine the geometry of the mag- 
netic field that the saturating non-linear dynamo gen- 
erates since this could be relevant for the infiuence of 
the field on e.g. asymmetric dust and wind formation. 
Qualitatively speaking the field becomes concentrated into 
elongated structures much thinner than the scale of the gi- 
ant convection cells, but perhaps due to the very irregular 
nature of the convective flows, no "intergranular network" 
is formed in the solar sense. On the one hand, at times 
magnetic structures coincide with downflows, but not as a 
general rule. On the other hand, strong fields are seldomly 
located within the general upflow regions. 

Figure |H1 shows the PDF of the magnetic field: the 
distribution is a typical signature of highly intermittent 
structures, i.e. only a very small fraction of the volume 
carries the strongest structures and the probability of find- 
ing a vanishing field strength at a random point in space, 
is far greater than finding strong fields. 



6 



Dorch: Magnetic activity in late-type giant stars 




Fig. 8. An illustration of the unsigned magnetic field strength |B| at the spherical surface r = R of the model star 
using an orthographic map projection. The darkest patches correspond to a maximum field strength of 500 Gauss 
(black on the continuous scale bar). From a snapshot at time = 695 years. The views are centered on longitudes of 
100° (left) and —100° (right). The grid indicated has a longitudinal spacing of 40° and a latitudinal spacing of 20°. 
The numerical resolution of the map is 180^ grid points. 




10 

Wavenumber (k) 

Fig. 7. Powerspectrum of the magnetic field B(fc)^ (solid 
curve) at time = 732 years, and a line corresponding to 
a power-law with an exponent of -2/3 (dashed line). The 
vertical lines, from left to right, denote the wavenumbers 
corresponding to the computational box size feg, the stellar 
radius fcrj,, and the Nyquist frequency fcNy (the numerical 
resolution) . 



An energy spectrum B^(fc) reveals that the mag- 
netic structures are well resolved with little power at the 
Nyquist wavenumber /cNy — tt/Ax, and that the power at 
the largest wavenumbers k ~ 100 is two-orders of magni- 
tude smaller than that at the largest scales (see Figure[7I). 
Maximum power is obtained on the largest scales corre- 
sponding to wavenumbers of a few, while there is a dip at 
fc ~ 7 corresponding to the scale of the radius, where the 




Fig. 9. A 3-d volume rendering of magnetic field lines 
(white) and flux ropes (dark structures). Also show is 
a isosurface (transparent) at a the surface temperature 
value. 



power is minimum. The power on scales k « 10-20 is flat 
leaning towards being proportional to fc^^/"^ correspond- 
ing to Kolmogorov scaling, and at small-scales fc ^ 50 
power steeply drops: magnetic structures in the non-linear 
regime are then large by solar standards, but smaller than 
the giant convection cells that show increasing power to- 
wards large scales. 
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Figure IHl is a map of the spherical surface at r = R in 
terms of magnetic energy in the non-hnear regime: there 
are both bright patches of strong magnetic field, e.g. a 
500 Gauss region at longitude around 80° and 50° north 
of the equator, and large areas with a vanishing field (e.g. 
at longitude 80° on the equator). On this map the areal 
magnetic filling factor is ~55% for |B| > 50 Gauss, while 
it is only ~0.6% for |B| > 500 Gauss. However, this map- 
ping does not represent any physical surface of the star. 
Due to the fact that the actual upper boundary consists of 
a few large cells the surface cannot be captured by a sim- 
ple sphere with radius R; this is illustrated by a volume 
rendering (Fig. (S)) of the isosurface at the cooling temper- 
ature. E.g. in the upper left corner of Figure there is 
a "hill" in the temperature isosurface and further large 
"slopes" can be seen across the star in this illustration. 
In the latter figure, the magnetic field lines illustrate that 
there is a slight trend inside the star, looking through the 
partly transparent surface, towards a radial orientation of 
the magnetic field, while the strong fields near the surface 
of the star are predominantly horizontally aligned. This 
was also observed by Dorch & Freytag (^2002^) and may 
be a generic trade of giant-cell full-radius convection in 
slowly rotating stars. 

4. Summary and conclusion 

In summary three different modes of dynamo action are 
recognized: 

1. A relatively fast growing linear mode with an expo- 
nential growth time of 4 years. 

2. A relatively slowly growing pseudo- linear mode with 
an exponential growth of ~ 25 years. 

3. A saturated non-linear mode operating a factor of two 
above equipartition. 

More modes may of course exist but these must then have 
very low growth rates and/or very small initial amplitudes 
since they have not appeared in the simulations. It is worth 
noting that in case 2) of the pseudo-linear mode, the same 
value of the growth time (around 25 years) was found in 
the previous purely kinematic dynamo models Dorch & 
Freytag (2002) although they employed a different com- 
putational method. This may in fact not be so strange, 
since the growth rate in a kinematic dynamo is set in part 
by the convergence of the flow across the field lines — -u 
and if the flows are similar so should the growth rates be. 

Based on the results presented here, it is not possible 
to state conclusively if Betelgeuse actually has a magnetic 
field, since such a field is unobserved. However, one may 
conclude that it seems possible that late-type giant stars 
such as Betelgeuse can indeed have presently undetected 
magnetic fields. These magnetic fields are likely to be close 
to or stronger than equipartition; they may be difficult 
to detect directly, due to the relatively small filling fac- 
tors of the strong fields, but even the moderately strong 
fields may have influence on their immediate surround- 
ings through altered dust, wind and mass-loss properties. 



The formation of dust in the presence of a magnetic field 
will be the subject of a subsequent paper along the lines 
presented here: the "Pencil Code" has recently been aug- 
mented with modules for radiation and dust modeling. 

The dynamos of the late-type giant studied here may 
be characterized as belonging to the class called "local 
small-scale dynamos" another example of which is the 
proposed dynamo action in the solar photosphere that is 
sometimes claimed to be responsible for the formation of 
small-scale flux tubes (cf. Cattaneo ll999f) . However, in the 
case of Betelgeuse this designation is less meaningful since 
the generated magnetic field is both global and large-scale, 
but because of the slow rotation, no large-scale solar-like 
toroidal field is formed. 

It is interesting to note that very recently Lobel et al. 
1)2004(1 published spatially resolved spectra of the upper 
chromosphere and dust envelope of Betelgeuse. Based on 
various emission lines they provide evidence for the pres- 
ence of warm chromospheric plasma away from the star 
at around 40 R. The spectra reveal that Betelgeuse's up- 
per chromosphere extends far beyond the circumstellar 
envelope. They compute that temperatures of the warm 
chromospheric gas exceed 2600 K. The presence of a hot 
chromosphere lead this author to speculate on the possible 
connection to coronal heating in the Sun, which is likely 
to be magnetic in origin and caused by flux braiding mo- 
tions in the solar photosphere (cf. Gudiksen & Nordlund 
I2()()2|l : it remains to be proven whether a similar process 
could be operating in late- type giant stars. 
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